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ABSTRACT 

We present the results of our analyses of the X-ray emission and of the strong lensing systems in the galaxy cluster Abell 611 
[z = 0.288]. This cluster is an optimal candidate for a comparison of the mass reconstructions obtained through X-ray and lensing 
techniques, because of its very relaxed dynamical appearance and its exceptional strong lensing system. We infer the X-ray mass es- 
timate deriving the density and temperature profile of the intra-cluster medium within the radius r^700 kpc through a non-parametric 
approach, taking advantage of the high spatial resolution of a Chandra observation. Assuming that the cluster is in hydrostatic equi- 
librium and adopting a matter density profile, we can recover the total mass distribution of Abell 611 via the X-ray data. Moreover, 
we derive the total projected mass in the central regions of Abell 611 performing a parametric analysis of its strong lensing features 
through the publicly available analysis software Lenstool. As a final step we compare the results obtained with both methods. We 
derive a good agreement between the X-ray and strong lensing total mass estimates in the central regions where the strong lensing 
constraints are present (i.e. within the radius r^lOO kpc), while a marginal disagreement is found between the two mass estimates 
when extrapolating the strong lensing results in the outer spatial range. We suggest that in this case the X-ray/strong lensing mass 
disagreement can be explained by an incorrect estimate of the relative contributions of the baryonic component and of the dark matter, 
caused by the intrinsic degeneracy between the different mass components in the strong lensing analysis. We discuss the effect of 
some possible systematic errors that influence both mass estimates. We find a slight dependence of the measurements of the X-ray 
temperatures (and therefore of the X-ray total masses) on the background adopted in the spectral analysis, with total deviations on the 
value of M 2 oo of the order of the lcr statistical error. The strong lensing mass results are instead sensitive to the parameterisation of the 
galactic halo mass in the central regions, in particular to the modelling of the Brightest Cluster Galaxy (BCG) baryonic component, 
which induces a significant scatter in the strong lensing mass results. 

Key words. Gravitational lensing: strong - X-ray: Galaxies: clusters - Galaxies: clusters: general - Galaxies: clusters: individual 
(A611). 



1. Introduction 

Testing the accuracy of the several approaches used to estimate 
the mass of galaxy clusters is necessary to assess the reliability 
of methods involving clusters as cosmological probes. For exam- 
ple, the comparison of the cluster baryon fractions fb to the cos- 
mic baryon fraction can provide a direct constraint on the mean 
mass density of the Un iverse, Qm (lEttori et al.ll2009l: I Allen et all 
120081: lEttori et al.ll2003h . while the evolution of the cluster mass 
function can tightly c onstrain Qa and the d ark energy equation- 
of-state parameter w ( Vikhl inin et al.|[2009b . 
Cluster mass profiles can be probed through several indepen- 
dent techniques, relying on different physical mechanisms and 
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requiring different assumptions. The comparison of results ob- 
tained with different methods and the (dis)agreement between 
them provide a useful check, and can give additional insights 
on the inner structure of clusters. For example, cluster masses 
can be estimated through the projected ph ase- space distribution 
and velocity dispersion of cluster galaxies (iBiviano & P oggianti 
20091; iDiaferio et al1 [2005): the dynamical studies, however, are 
affected by the strong degeneracy between the mass density pro- 
file and the velocity dispersion anisotropy profile, and thus re- 
quire strong assumptions on the form of the velocity dispersion 
tensor. 

Measurements of the Sunyae v-Zerdovich effect (hereafter SZE; 
ISunyaev & Zerdovichll 19721) towards galaxy clusters yield a di- 
rect measurement of the pressure distribution of the cluster gas: 
the total gravitational mass can thus be determined combining 
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Fig. 1 : [Left panel] F606W ACS/HST image of the cluster Abell 611 (see Table[3]for the observation summary). The X-ray isophotes 
derived from the Chandra image (shown in the middle panel) are overlaid in white. The green cross indicates the X-ray centroid. 
[Middle panel] The Chandra X-ray image of Abell 611 in the energy range [0.25-7.0] keV. The image was smoothed with a 3 
pixel FWHM Gaussian filter, and is not exposure-corrected. The contours derived from the HST image (shown in the left panel) 
are overlaid in white to facilitate the comparison between the images. Here again the green cross marks the X-ray centroid. [Right 
panel] Pseudo-colour composite image of Abell 611, obtained combining r and g-SLOAN images, taken with the Large Binocular 
Camera mounted at LBT. The three panels are WCS-aligned; the size of the field of view is ^ [1.6 x 1.8] arcmin. 



this information with the intra-cluster medium (hereafter ICM) 
temperatur e, estimated for e xample through the X-ray spectral 
analysis (iGrego et al.l l200ll) . Because it is independent of the 
cluster redshift, the SZE can be an extremel y powerful tool to in - 
vestigate the mass of high-redshift clusters (Hi ncks et al.ll2009l) . 
but accurate measurements of the SZE are still a challenge be- 
cause of intrinsic limitations and systematic biases (faint signal 
extending over a large angular size and subject to several sources 
of contamination). 

So far, the analysis of the cluster X-ray emission and of the grav- 
itational lensing effect are among the most promising techniques 
to estimate galaxy cluster masses. 

The X-ray mass estimate is less biased compared with lensing- 
derived masses by projection effects, because the X-ray emission 
is tra ced by the square of the gas density (iGavazzil [20051: IWul 
2000). However, X-ray measurements of cluster masses imply 
the assumption of hydrostatic equilibrium of the ICM with the 
dark matter potential and of spherical symmetry of the cluster 
mass distribution. Hence the total mass profile can be inferred 
from the radial profiles of gas temperature and density. 
On the other hand, the gravitational lensing effect allows for the 
determination of the projected surface mass density of the lens, 
regardless of its dynamical state or the nature of the intervening 
matter. This effect is determined by all massive structures along 
the line of sight, so lensing mass measurements are subject to 
foreground and background contaminations. Moreover, the lens- 
ing effect is sensitive to features of the mass distribution such as 
its ellipticity and asymme try and to the presence of substructures 
(Me neghetti et al.ll2007ah . 

The combination of X-ray and lensing results could verify and 
strengthen th e findings on the cluster physics and masses (see 
for example lAUenlll998b lEttori & Lombardill2003l: iBradac et al.l 
l2008h : but so far several works in literature claimed a signifi- 
cant disagr eement between strong lensing and X-ray mass es- 
timates, (IWu & Fand |T99g; ISmail et all 119971: lota et al I [2004; 
IVoigt & Fabianl 120061: iGitti et al.l 120071: lHalkola et al.l [2008) 
which would prevent any joint analysis. Many convincing ex- 



planations, even if not conclusive, have been suggested. 
It is therefore clear that a comparison between independent mass 
estimates for galaxy clusters can provide a vital check of the as- 
sumptions adopted in the cluster analysis, and will possibly give 
us additional insights into the underlying physics of these ob- 
jects. In particular, the comparison between lensing and X-ray 
mass in galaxy clusters is also directly related to the cosmo- 
logical parameters, since their ratio depends on the cosmic den- 
sity parameter, Qm, and the cosmological constant, Qa, for their 
different dependence on the angular diameter distances, w hich 
are s maller in an Qa dominated universe, as pointed out by IWul 
d2000b . 

Targets for a comparative analysis should be relaxed, unper- 
turbed objects to avoid biases deriving from the incorrect as- 
sumption of hydrostatic equilibrium (X-ray analysis) or from 
the contamination of secondary mass clumps or unresolved sub- 
structures (lensing analysis). 

An ideal candidate for this kind of an alysis is Abell 611. It is a 
rich, cool-core cluster at z = 0.288 dStruble & Roodlll999l) : its 
X-ray emission peak is well centred on the Brightest Cluster 
Galaxy (hereafter BCG; see left and central panels of Fig. [U, 
and the X-ray isophotes appear quite regular, with a low de- 
gree of ellipticity. These characteristics are generally considered 
good indicators of a relaxed dynamical state. Relaxed objects are 
quite rare at this redshift, so a detailed study of the properties of 
Abell 61 1 is even more appealing. 

In this work we will focus on the comparison between the X-ray 
and strong lensing mass estimates. A weak lensing analysis 
of large field of view data, obtained with the Large Binocular 
Camera (LBC) mounte d at the Lar g e Bin ocular Telescope 
(LBT), is presented in iRomano et al.l (1201 Oh (hereafter Paper 
I). Here we present new X-ray and strong lensing mass esti- 
mates, based on a non-parametric reconstruction of gas den- 
sity and temperature profiles obtained through the analysis of 
a Chandra observation, and on a strong lensing reconstruction 
respectivel y, which was performed with the Lenstool analysis 
software (Kne ib et alJfl995 lJullo et alJl2007h . 
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Table 1 : Observation summary of the Chandra exposure used for 
the X-ray analysis. 



Obs.ID 


Start Time 


Total 


Net 


PI 






Expos, [ks] 


Expos, [ks] 


name 


3194 


2001-11-03 18:43:58 


36.6 


30.9 


Allen 



This paper is organised as follows. In Sect. [2] we will intro- 
duce our X-ray analysis, focusing on the data reduction and 
on the method applied to recover the total and gas mass pro- 
files, and discuss some possible systematic errors associated 
with the X-ray analysis. The strong lensing analysis is presented 
in Sect. where we will briefly discuss our main findings. We 
will compare the X-ray and strong lensing results in Sect. HI 
a comparison with previous analyses can be found in Sect. 
Finally, we will summarise our results and draw our conclusions 
in Sect. [6] 

Throughout this work we will assume a flat ACDM cosmol- 
ogy: the matter density parameter will have the value £_ m =0.3, 
the cosmological constant density parameter Qa=0.7, and the 
Hubble constant will be Ho = 70 km s _1 Mpc _1 . At the cluster 
redshift and for the assumed cosmological parameters, 1 arcsec 
is equivalent to 4.329 kpc. Unless otherwise stated, all uncertain- 
ties are referred to a 68% confidence level. 

2. X-ray analysis 

2.1. X-ray data reduction 

The X-ray analysis of Abell 611 returns several indicators of 
a relaxed dynamical state, such as the absence of substruc- 
tures in the X-ray emission (which appears isotropic over ellip- 
tical isophotes), a well-defined cool-core, and a smooth surface 
brightness profile. The X-ray emission is centred on the BCG: 
the distance between the X-ray centroid and the BCG centre is 
=- 1 arcsec (the uncertainty on this measure is comparable to 
the smoothing scale applied to the X-ray image to determine the 
centroid, i.e. 3 pixels ^1.5 arcsec). Moreover, Abell 611 ap- 
pears to be a radio-quiet cluster, because the Gia nt Metrewave 
Radio Teles cope (GMRT) Radio Halo Survey (iVenturi et al.l 
120071 120081 see Sect. 0) failed to detect any extended radio 
emission associated with this cluster. Since giant radio halos are 
the most relevant examples of non-thermal activity in galaxy 
clusters (iBrunetti et al.l 12009 ), the absence of a radio halo sup- 
ports the hypothesis of thermal pressure support, connected to 
the hydrostatic equilibrium assumption required in the X-ray 
analysis. 

This cluster can therefore be considered an ideal candidate to 
highlight any issue with the current analysis techniques: it full- 
fills the assumptions of thermal pressure support and of spher- 
ical symmetry, so it is unlikely that potential biases in the re- 
sults are related to the limiting assumptions underlying the X-ray 
analysis. We performed our X-ray analysis on the only Chandra 
data set available, retrieved from the public archive (see Table Q] 
for observation log). We reduced the data using the CI AO data 
analysis package - version 4.0 - and the calibration database 
CALDB 3.5.(fJ We will summarise here briefly the reduction 

1 See the CIAO analysis guides for the data reduction: 
cxc.harvard.edu/ciao/guides/ . 



procedure. 

The Chandra observation was tele-metered in very faint mode 
and was performed through both back-illuminated (BI) and 
front-illuminated (FI) chips. The cluster centre was imaged 
in the S3 BI chip. The level- 1 event files were reprocessed 
to apply the appropriate gain maps and calibration products 
and to reduce the ACIS quiescent background We used the 
acis_process_events tool to check for cosmic-ray back- 
ground events and to correct for eventual spatial gain varia- 
tions caused by charge transfer inefficiency to re-compute the 
event grades. We applied the standard event selection, including 
only the events flagged with grades 0, 2, 3, 4, and filtered for the 
Good Time Intervals associated with the observation. The bright 
point sources were masked out, and a geometrical correction for 
the masked areas was applied; the point sources were identified 
using the script vtpdetect, and the result was then checked 
through visual inspection. 

A careful screening of the background ligh t curve is neces- 
sary for a correct background subtraction (iMarkevitch et al.l 
120031) . to discard contaminating flare events, which could have 
a non-negligible influence on the inferred cluster spectral prop- 
erties. We extracted the background light curve, applying a time 
binning of about 1 ks, in the energy range [2.5-7] keV, which is 
the most effective band to individuate common flares for S3 chip. 
We applied the script lc .clean to include only the time periods 
within a factor 1.2 of the mean quiescent rate. We compared the 
S3 background light curve with the light curve extracted in the 
SI chip using the energy range [2.5-6] keV, to check for faint, 
long flares, which were not identified. The S3 background light 
curve was examined using the ChlPS facilities to identify and 
exclude further flaring events. The net exposure time after the 
light curve screening is about 31 ks (see Table [T]). 

2.2. X-ray analysis 

In order to derive a reliable cluster mass estimate through the 
X-ray analysis, a primary issue is to avoid any bias in the tem- 
perature estimate, because cluster masses derived by assuming 
hydrostatic equilibrium are dependent on the temperature pro- 
file, as demonstrated by iRasia et al.l (120061) . For this purpose, a 
correct background modelling is crucial. The X-ray background 
is often estimated through the blank- sky background data sets 
provided by the ACIS calibration team. The blank-sky observa- 
tions are reduced similarly to the source data sets to be compati- 
ble with the cluster event files; the former are then re-normalised 
by comparing the blank- sky count rate in a given energy range 
(generally kT > 8.0 keV, since in this band the Chandra effec- 
tive area is negligible) with the local background count rate. One 
of the advantages of this modus operandi is that the derived ARF 
and RESPONSE matrices will be consistent both for the source 
and the background spectrum. However, the background in the 
X-ray soft band can vary both in time and in space, so it is impor- 
tant to verify whether the background derived by the blank- sky 
datasets is consistent with the "real" one. 

For this purpose, we extracted a spectrum from the Chandra ob- 
servation in a region where the cluster emission is negligible: 
we compared the former to a spectrum derived in the same re- 
gion of the blank-field data set. We tried to fit the residuals in 
the [0.4-3] keV band with a MEKAL model, without an absorp- 
tion component and broadening the normalisation fitting range 
to negative values. This comparison shows some differences be- 



2 For a complete discussion on this topic, see 
cxc.harvard.edu/cal/Acis/Cal_prods/vfbkgrnd/index.html . 
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tween the two spectra. For this reason, we decided to proceed in 
the spectral fitting utilising both a local background estimate and 
a blank- sky estimate. The effect of the background choice in the 
mass estimate results will be discussed in Sect. 12.51 



2.3. Spectral fitting 

The spectra were extracted on concentric annuli centred on the 
X-ray surface brightness centroid, each one containing at least 
2500 net source counts, to infer an estimate of the ICM metallic- 
ity. We selected six annuli with adequate S/N ratio, up to R spQC = 
2.9 arcmin. We used the CIAO specextract tool to extract 
the source and background spectra and to construct ancillary- 
response and response matrices. The spectra were fitted in the 
0.6 - 7.0 keV range except for the last two annuli, in which, 
due to the higher background level, we restricted the analysis to 
the range 0.6 - 5.0 keV. We perfor med the spectr al analysis us- 
ing the XSPEC software package (lArnaudlll996l) . We adopted 
an optically-thin plasma emission model (the mekal model; 
iMewe et all 1 19851: iKaastra & Mewd 1 19931) with a n absorption 
comp onent (the Tubingen-Boulder model tbabs; Iwilms et all 
l2000h . 

We fixed t he Galactic absorption to the value inferred from radio 
HI maps in lDickev & Lockmanl(ll990h . i.e. 4.99xl0 20 cm" 2 . The 
free parameters in the spectral fitting model are the temperature, 
the metallicity and the normalisation of the thermal spectrum. 
The metal abundance profile presents a central peak with a quite 
scattered trend (see the right panel of Fig. [2]). The metallicity 
profile is not directly involved in the X-ray mass reconstruction: 
but we still verified whether a bias in its measurement can influ- 
ence the derived mass estimate. The results of this check will be 
presented in Sect. 



2.4. X-ray mass profile 

The gas temperature and density profiles are recovered without 
any parameterisation of the ICM properties, but relying on the 
assumptions of spherical symmetry and hydrostatic equilibrium 
and on the choice of a model for the total mass. 
The method that was app lied is briefly resumed in the follow - 
ing; the reader can refer tolEttori. De Grandi, & Molendi ( 2002); 
Mor andi. Ettori. & Moscardinil(l2007h :lMorandi & Ettoril d2007) 
for full details. 

In this section, the subscript "ring" will indicate a projected 
quantity measured in concentric rings (equivalently, spectral 
bins), whereas the subscript "shell" will point to a deprojected 
quantity measured in volume shells. The projected gas temper- 
ature r r i n g, metal abundance Z r i ng and the normalisation K Y i ng of 
the model for each annulus are derived through the spectral anal- 
ysis. In particular, the normalisation of the thermal spectrum pro- 
vides an estimate for the emission integral EI = J n e n p dV ^ 
0.82 Jrc 2 dV. 

To deproject the former quantities, we need to derive the inter- 
section of each volume shell with the spectral annuli, which can 
be calculated through an upper triangular matrix r V, whose last 
pivot represents the outermost annulus if we assume spherical 
symmetry dKriss. Cioffi. & Canizares1ll983UBuotell200of) . 
To calculate the deprojected density profile, we combined the 
spectral information and the cluster X-ray surface brightness 
(SB) profile, composed by 15 bins with 1000 counts each. The 
SB provides an estimate of the volume-counts emissivity oc 
n ebin^b/n /^L' wnere A. is the luminosity distance of the cluster. 
By comparing this observed profile with the values predicted by 



a thermal plasma model, with temperature and metallicity equal 
to the measured spectral quantities and taking into account the 
absorption effect, we can solve for the electron density and ob- 
tain a gas density profile better resolved than a "spectral-only" 
one (see the top panel of Fig. [2]). 

The total gravitating mass was then constrained as follows: the 
observed deprojected temperature profile in volume shells r s h e ii 
was obtained through 



6r shell = (T r )- 1 #(L ring r ring ), 



where 



e = (T ) #Lj 



Ting > 



(i) 



(2) 



here the symbol # indicates a matrix product, *V is the volume 
matrix and L r i ng is the cluster luminosity measured in concentric 
rings. We compare it to the predicted temperature profile r mo deb 
obtained by inverting the equation of hydrostatic equilibrium be- 
tween the dark matter and the intracluster plasma, i.e.: 



- Gjim ] 



fleMot,model(< r) _ d (ft e X ^model) 

r 2 dr 



(3) 



where n Q is the deprojected electron density. Our best- fit mass 
model was determined by comparing the observed temperature 
profile r bs with the predicted profile r mo dei> which depends on 
the mass model parameters: accordingly we minimised the quan- 
tity: 



X 



rings 



model 



0~ obs 



(4) 



where the r mo d e i profile was rebinned to match the observed 
temperature profile. In the bottom panel of Fig. [2] we show the 
projected temperature profile as determined through the spectral 
analysis, the deprojected best-fit profile and the deprojected pro- 
file, rebinned to the spectral intervals. 

We assumed the NFW profile as the total mass functional; this 
profile can be expressed as: 



P(r) = 



Ps 



(r/r s )(l+r/r 5 Y 



(5) 



where p s is the characteristic density and r s is the scale radius. 
We will describe the NFW profile through the scale radius and 
the concentration parameter C200; me former can be defined as 



R 



c = 



200 



(6) 



where R200 is the radius within which the mean matter density 
of the halo is 200 times the critical density of the Universe. 
Assuming a mass model, the best-fitting values of its parame- 
ters are determined by minimising the quantity defined in Eq. [4J 
The NFW profile parameters were optimised in the ranges 0.5 < 
C200 ^ 20 and 10 kpc < r s < 976 kpc, where the scale radius up- 
per limit is the outer radius of the surface brightness profile. The 
gas and total mass profiles are shown in Fig.[3l the inferred mass 
model parameters are listed in Table 12 The confidence levels for 
the NFW parameters are shown in the right panel of Fig. [3] 

2.5. Systematic effects in the X-ray mass profile 

In this section we will briefly assess if our X-ray mass estimate 
of Abell 611 can be affected by some analysis choices. 
The choice of the X-ray background estimate, as discussed in 
Sect. 12.21 could have a non-negligible impact on the derived 
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Fig. 2: [Top left panel] Deprojected gas density profile. The red circles indicate the values determined by the normalisation of the 
thermal spectra, while diamonds represent the values obtained combining both spectral and surface brightness measurements: the 
dashed line shows the best-fit profile. [Bottom left panel] ICM temperature profile. The figure shows the projected temperature 
values, directly inferred from the spectral analysis (red filled circles) and the deprojected values, which were rebinned to the same 
binning of the spectral profile (diamonds). The dashed line shows the deprojected best-fit profile. [Right panel] The projected (black 
empty diamonds) and deprojected (red filled circles) metallicity profile for Abell 611 inferred from the X-ray analysis. The dotted 
line represents the mean metallicity value Z = (0.45 ± 0.1) Z Q , as determined through the fit of the spectrum extracted in the area 
with radius R= 2.9 arcmin, i.e. the cluster area covered by our spectral analysis. The grey region indicates the l<x confidence level. 

Table 2: Best-fit values for the total mass parameters obtained through the X-ray analysis. 



Results of the X-ray analysis of Abell 611 





r s 


C200 


M200 




o.f.] 




[kpc] 




[1O 14 M ] 






(1) Blank- sky background 


350.3 + 79.6 


5.18 + 0.84 


9.32+1.39 


0.81 


[4] 


(2) Local background 


420.0 + 99.5 


4.67 + 0.80 


11.11 ±2.06 


0.46 


[4] 


(3) Fixed metallicity 


367.2 + 98.6 


5.13 + 0.87 


10.61 + 1.96 


0.33 


[4] 



Notes. (1) the blank- sky dataset is adopted to estimate the X-ray background; 

(2) the X-ray background is derived from peripheral regions of the source observation; 

(3) the background is estimated locally fixing the gas metallicity to Z = 0.3 Z©. 

The columns indicate the best-fit values obtained for the NFW profile parameters (the scale radius r s and the concentration C200X me estimate of 
the total mass within 7? 2 oo (here listed as M 2 oo) and the reduced^ 2 , referred to a fit with four degrees of freedom. 



temperature profile, and consequently on the total mass mea- 
surement. Moreover, we verified whether considering the ICM 
metallicity as a free parameter in the spectral analysis can im- 
pact the X-ray mass estimate of Abell 611, since (as shown in 
Fig. [2]) its inferred metallicity profile is quite irregular. We thus 
re-ran the spectral analysis by fixing the met allicity Z to the typ- 
ical value of 0.3Z© (see e.g. Bal estra et al. 112007b : we then de- 
rived a new X-ray mass estimate adopting the spectral profiles 
obtained under the assumption of a constant metallicity. 
Therefore we derived the mass of Abell 611 through the X-ray 
analysis: 



1 . by estimating the X-ray background through the spectra ex- 
tracted from the blank-sky dataset (opportunely reprocessed 
and re-normalised to be consistent with the source observa- 
tions) in the same image regions used for the source spectra 
(case 1); 

2. by estimating the X-ray background from peripheral regions 
of the S3 chip in the source dataset (case 2); 

3. as in the previous case, but by fixing the metallicity Z to the 
value 0.3Z o (case 3). 

The NFW parameters derived from the X-ray analysis in the 

three aforementioned cases are listed in Table [2 

The results obtained in the first two cases mutually agree within 
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Fig. 3: [Left panel] Mass profiles for Abell 611 as determined through the X-ray analysis. The black solid line and the red dashed 
line indicate the total mass and the gas mass profile, respectively. [Right panel] Confidence levels determined through the X-ray 
analysis for the NFW profile parameters r s and C200. The red circle indicates the best-fit values, while the contours represent the 
lcr, 2<t, and 3<x confidence levels. 
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Fig. 4: Comparison between the temperature profile determined 
by estimating the X-ray background through the blank- sky 
dataset (red dashed line) or using a local background component 
(black solid line). 



the lcr range, even if the temperature profile derived in the latter 
case is systematically higher (see Fig. |4j). Moreover, the best-fit 
parameters obtained in the third case (i.e. fixing the metallicity 
to 0.3Z o : see case 3 in Table [2]) agree with the values obtained 
by including the metallicity as a free parameter in the fitting: the 
choice of estimating or fixing the metallicity in the spectral anal- 
ysis has a negligible impact on the derived cluster mass. 



3. Strong lensing analysis 

3.1. Imaging data 

Our strong lensing analysis of Abell 611 benefits from the 
high-resolution imaging data obtained with the Hubble Space 



Table 3: Observation summary of the ACS/HST image used for 
the strong lensing analysis. 



Data-set ID 
Start Time 
Exposure [ks] 
Instrument 
Filter 

Proposal ID 
PI 



J8D102010 
002-12-03 21:07:28 
2.16 

ACS/WFC 
F606W 
9270 
Allen 



Telescope (HST), which allowed an unambiguous identifica- 
tion of several lensed systems. Abell 611 was observed with 
the Advanced Camera for Surveys (ACS) mounted on HST. 
We retrieved the ACS association data of Abell 611 from the 
Multimission Archive at STScI (MASTfl; the details of the HST 
observation of Abell 61 1 are summarised in Table [3j 
The HST archival data for Abell 611 are already processed fol- 
lowing th e ACS pipeline and com bined using the multi drizzle 
software dKoekemoer et al.l2002l) . with a final pixel scale of 0.05 
arcsec; consequently no additional basic data reduction is neces- 
sary. But we used the soft wares SExtra ctor dBertin & Arnoutsl 
119961) . Scamp (lBertinll2006h and Swar$ dBertin et alJ2002h . de- 
veloped at the TERAPIX data reduction centre to derive the ob- 
ject catalogue for the Abell 611 HST imaging data, re-compute 
an astrometric solution and finally resample the image to de- 
rive a more accurate astrometric matching and to equalise the 
background level. The WCS coordinates reported in the present 
work are given with regard to the HST data and refer to the as- 
trometric solution that we re-derived. Abell 611 was also ob- 
served during the Science Demonstration Time for the Blue 
Channel of the Large Binocular Camera (LBC), mounted at the 
Large Binocular Telescope (LBT) (located at the Mt. Graham 



3 The url of the MAST is archive.stsci.edu. 

4 The software SExtractor, Scamp and Swarp can be found at the url: 
astromatic.iap.fr/software/sextractor, astromatic.iap.fr/software/scamp 
and astromatic.iap.fr/software/swarp, respectively. 
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International Observatory, Arizona), in good seeing conditions 
(FWHM of the seeing disc ~ 0.6 arcsec). The observations 
were made with the g-SLOAN and r filters; the total exposure 
times of the two coadded images are 3.6 ks and 0.9 ks, respec- 
tively. The LBT data were reduced according to the data re- 
duction pipeline developed at the INAF/OAR centre^; the as- 
trometric solution wa s derived with the ASTROMC software 
dRadovich et al.ll2008h . The final plate scale of the LBC images 
is 0.225 arcsec/pixel. Further details on the LBC/ LBT data and 
on the ir reduction can be found in Paper I and in Giallo ngo et al.l 
(2008); a comparison with the weak lensing mass estimate ob- 
tained in Paper I is presented in Sect. The pseudo-colour im- 
age of Abell 611, obtained by combining the two LBT observa- 
tions, is shown in Fig.[TJ We relied on the g-SLOAN and r LBT 
observations both to verify whether the candidate conjugated im- 
ages have similar colours, as it should be because of the intrinsic 
surface brightness conservation in gravitational lensing, and to 
select the cluster galaxies , in order to account for their mass in 
the lens modelling (see § 13.51) . 




3.2. Strong lensing system 

Abell 611 exhibits an outstanding strong lensing system. The 
most evident lensed feature is the giant tangential arc located at 
about 15.5 arcsec from the BCG centre (see Fig. [5]). The cur- 
vature of the arc significantly changes in the proximity of some 
cluster galaxies, which locally perturb the critical lines and pos- 
sibly t he magnification and the multiplicity of the lensed im- 
ages dKee ton 2003 llAmara et al . 2006; Mene ghetti et alJl2007bt 
Peirani et al. 2008; Puchwein & H ilbert 20Q9|) 
Richard et al I d2010h (hereafter R09) obtained the optical spectra 
of three surface brightness spots along the tangential arc, deriv- 
ing a redshift estimate of z = 0.908 ± 0.005 for all of them. The 
arc appears to be formed by three merging images (A.l, A.2 and 
A.3 in Fig. O, as in a typical cusp configuration. 
Another confirmed lensed system is a fold system (hereafter 
called system B). iRichard et al.l (l201Qh reported the spectro- 
scopic redshifts obtained for images B.2 and B.3 in Fig. [3 con- 
firming that both of them are lensed images of a source at red- 
shift z = 2.06 ± 0.02. The images labelled B.l and B.4/B.5 can 
be safely considered as their counter-images. More specifically, 
the image B.5 appears to be an additional image produced by the 
position of the source B, lying totally or partially inside the caus- 
tic of the nearby galaxy (galaxy n.l in Table @]). The g-r colours 
of all images of system B are similar (see Fig.Q]). 
A third lensed system was serendipitously identified by R09. 
They were able to derive the redshift zc.i = 2.59 ± 0.01. of im- 
age C.l (see Fig.0 thanks to its strong Lyman-alpha emission. 
The image C.2 (see Fig. was considered by R09 to be the 
counter-image of C.l. We agree with their identification and in- 
cluded both images as constraints. 

Finally, the images labelled D.l to D.4 in Fig. \5\ were identi- 
fied by R09 and bv iNewman et al.l (l2009h (hereafter N09; see 
Sect. [5] for more details about their recent study of Abell 611) 
as a lensed system, from the similarity of the morphologies of 
images D.l and D.2. However, this system is not yet confirmed 
spectroscopically, so we decided not to include it in the strong 
lensing analysis, but we verified if and how the best-fitting re- 
sults would be modified by considering it as an additional con- 
straint (see Sect. 13.71) . 



5 See the Large Binocular Camera website hosted by the OAR- 
Osservatorio Astronomico di Roma: lbc.mporzio.astro.it . 



Fig. 5: Overview of the strong lensing system of Abell 611. The 
knots that we used as constraints are marked as red crosses on the 
F606W ACS/HST image; enlarged images are shown in Fig. IA.ll 
The magenta crosses indicate the perturber galaxies that were 
optimised singularly (see Case 6 in Table [5]). The blue crosses 
correspond to a likely lensed system, identified by R09, which 
is not yet confirmed spectroscopically. We discuss our selec- 
tion of the input images and its effect on the results in Sect. 13.21 
and Sect. 13.71 The coordinates of all the systems are listed in 
Table IA. ll The size of the field of view is ^ 50 x 42.5 arcsec; 
the image was smoothed with a 2 pixel FWHM Gaussian filter 
(1 pixel = 0.05 arcsec). 



3.3. Strong lensing modelling 

We performed a parametric reconstruction of the mass dis- 
tribution of Abell 611 with the L enstool analysis softwar^ 
(iKneib et al.lll9%l:lJullo et al.ll2007h . 

Our lens model consists of the superposition of a smooth 
cluster-scale potential and of galaxy-scale halos, associated with 
early-type cluster galaxies. The d etails of the cluster scale halo 
modelling will be discussed in § 13.41 whereas we will discuss 
the galaxy- scale component in § 13.51 finally the BCG modelling 
will be discussed separately in § 13.61 

We constrained the lens model for Abell 611 by identifying in 
each image system, consisting of the multiple images of the 
same source, some elliptical regions or knots, which are likely 
to be multiple images of the same source area. The knots were 
identified through visual inspection of the ACS/HST observation 
by selecting the local peaks in the surface brightness, which are 
detectable in the strong lensing images of Abell 611. The surface 
brightness spots were then grouped into sets of conjugated knots 
by evaluating characteristics such as their relative positions with 
respect to the location of the critical lines and their brightness 
ratios. 

We assumed an astrometric positional error of <x p0 s = 0.2 arc- 
sec. The adopted error affects the^ 2 estimator computed in the 



6 The Lenstool software package is available at 
www. oamp . f r/co smology/lenstool . 
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source plane by Lenstool: 

where x J s is the position of the source corresponding to the ob- 
served image j, < x J s > the barycenter of the ni source positions 
(corresponding to the ni conjugated images or knots of the image 
system i), and fij the magnification for the lensed image j; the 
final x 2 is then obtained by combining the individual x\ com- 
puted for each image system. 

Between the several optimisation methods available in Lenstool 
we adopted the Bayesian optimisation method, with a con- 
vergence speed parameter equal to 0.4; this method, although 
slower than the other optimisation options, is less sensitive to 
local minima in the parameter space. 

3.4. Cluster-scale component 

We modelled the cluster scale component assuming an elliptical 
NFW profile as mass functional. The initial values for the 
centre, ellipticity e and position angle of the cluster halo 
were inferred from the BCG luminosity profile. The centre 
was allowed to vary ±5 arsec, while the ellipticity and the 
position angle were optimised in the ranges 0.01 < e < 0.4 and 
100.0 deg < #dm < 140.0 deg, respectively. The optimisation 
range of the cluster position angle is asymmetric in order to 
take into account the orientation of the X-ray isophotes. The 
optimisation ranges for the NFW profile parameters were 
C200 € [0.5 - 14.0] and r s e [50.0 - 800.0] kpc. 



3.5. Galaxy-scale components 
3.5.1 . Selection of cluster galaxies 

Several authors demonstrated that the inclusion of galaxy-scale 
halos in the lens modelling (in particular the galaxies whose 
projected position is within the region probed by strong lensing) 
is necessary to account for their effect on the lensing cross 
section and on the appear ance of the lensed images that fall 
close to the perturbe rs (Nataraia n et al.l I2007I; lAmara et~aD 
l2006UMeneghetti et al . 2003; Keeton 2003; Bradac et al. 2QQ3). 
Although gravitational lensing is sensitive to all massive struc- 
tures distributed along the line of sight, the main galactic mass 
budget is constituted of the massive early-type cluster galaxies; 
we will thus consider only this galactic population in our lens 
modelling. 

The cluster members were selected through the characteristic 
cluster red-sequence, which a l lows us to identify the red 
early-type galaxies (iBell et al.l [2004) and which is clearly 
detected in the colour-magnitude diagram (see Fig. 0. We 
will assume in the modelling that all the red-sequence galaxies 
belong to the cluster, and thus lie at its same redshift. The 
diagram was derived by extracting the object catalogues from 
the two LBT observations, using the SExtractor software in 
single and dual-image mode. 

Even if the total exposure time of the red image is 0.9 ks, which 
is one fourth of the g-Sloan observation exposure time, the 
former filter allows an easier deblending of the redder cluster 
galaxies from the bluer lensed features. For this reason we used 
the r filter image for the detection of the sources in the dual- 
image mode. Moreover, even if the higher spatial resolution 




mag R 



Fig. 6: Colour-magnitude diagram of Abell 611. The black lines 
define the cluster red-sequence, while the blue filled squares in- 
dicate the red sequence galaxies lying within 100 arcsec from 
the BCG centre, which were included in the lens modelling. The 
black circle indicates the galaxy located between the lensed im- 
ages B.l and C.l (see Fig. 0, which was included in the lens 
modelling for its position; the magnitude estimate derived for 
this galaxy from the LBT data is contaminated by the flux con- 
tribution of the nearby lensed feature (see text for details). 

of ACS data would guarantee a more accurate deblending of 
contiguous sources, we decided to use the LBC data, which 
cover a wider area, to derive the cluster galaxy catalogue. 
Only galaxies located within 100 arcsec (^ 433 kpc) from the 
BCG centre and whose magnitude Mr is lower than 23 were in- 
cluded in the lens modelling; a total of 97 galaxies were selected 
adopting these criteria. We added to the catalogue the galaxy 
located between system B.2 and B.3 in Fig. [5] that could affect 
the lensed images due to its position. The colour-magnitude 
diagram, the cluster red- sequence and the galaxies included in 
the lens modelling are shown in Fig. [6] 



3.5.2. Estimate of cluster galaxy parameters 

We used the tw o-dimensional fitting algorithm GALFIT^ 
dPeng et al.l l2QQ2h to model and subtract the BCG as well as 
seven galaxies located within ^17 arcsec from the BCG cen- 
tre and with magnitude in the ACS F606W filter lower than 22 
for a better identification of the multiple lensed images and to 
extract the BCG surface brightness profile. All galaxies were fit- 
ted simultaneously; the lensed features and the point-like sources 
were masked out. We modelled the galaxy light profile adopting 
both a Sersic and a de Vaucouleurs profile. The Sersic profile 
(ISersidll968l) has the functional form 









-K 


((;)""-) 











where H e is the pixel surface brightness at the effective radius 
r e , the latter defined as the radius that encloses half of the to- 
tal flux, and n is the Sersic index. The de Vaucouleurs profile 
(Ide Vaucouleurs! 1 948b has the same analytical form as the Sersic 

7 The GALFIT software is publicly available at the url 
user s . obs . carnegiescience . edu/peng/work/galfit/galfit. html . 
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Table 4: Best-fit parameters obtained by modelling the galaxy luminosity with the Sersic profile (ISersid [l968). 



Best-fitting Sersic profile parameters 



ID 










JT VJ5>111VJ11 


Sersic 


m 606W 




deg(J2000) 


deg(J2000) 


ratio 


[kpc] 


Angle 


index 




BCG 


120.23678 


36.056572 


0.70 + 0.01 


38.54 + 0.27 


132.5 + 0.1 


3.0 + 0.1 


17.0 + 0.1 


1 


120.23598 


36.054391 


0.83 + 0.01 


1.48 + 0.01 


112.8 + 0.1 


2.3 + 0.1 


20.7 + 0.1 


2 


120.2357 


36.053776 


0.92 + 0.01 


1.02 + 0.01 


13.1 + 0.1 


2.9 + 0.1 


21.6 + 0.1 


3 


120.23855 


36.061407 


0.50 + 0.01 


3.85 + 0.04 


78.7 + 0.2 


2.6 + 0.1 


20.7 + 0.1 


4 


120.24051 


36.059402 


0.67 + 0.01 


3.28 + 0.07 


80.6 + 0.5 


4.7 + 0.1 


20.8 + 0.1 


5 


120.24255 


36.056744 


0.84 + 0.01 


1.35 + 0.02 


128.4 + 1.4 


4.6 + 0.1 


21.4 + 0.1 


6 


120.23639 


36.055754 


0.90 + 0.01 


0.64 + 0.01 


131.6 + 4.0 


2.3 + 0.1 


22.0 + 0.1 


7 


120.24148 


36.060146 


0.79 + 0.01 


2.32 + 0.04 


61.7 + 0.9 


5.1+0.1 


21.1+0.1 



Notes. The algorithm GALFIT is used. The fit was performed using the ACS F606W filter image. We list here the galaxy ID (the galaxy numbering 
is the same used in Fig. [5}, the centroid coordinates, the axis ratio, the effective radius r e , the positional angle, the Sersic index n and the galaxy 
integrated magnitude, computed in the ST magnitude system. Here we report the 68% statistical errors as determined from GALFIT. The errors 
were rounded to the closest value greater than zero. 



profile, but the index has the fixed value of 4. 
The galaxy models were convolved with the HS T PSF image cre- 
ated using the Tiny Tim software (Kriii l 1993b . The reduced^ 2 
of the best-fitting models obtained adopting either the Sersic or 
the de Vaucouleurs profile are similar,^ =3.8 and^ 2 y = 3.9, 
respectively. However, the galaxy- subtracted image derived fit- 
ting the galaxy light profile with a de Vaucouleurs model shows 
larger residuals in the central area for some of the fitted galaxies, 
in particular for the galaxy labelled 1 in Table |4] For this reason 
we chose a Sersic profile to model and subtract from the HST im- 
age the galaxies located within ^17 arcsec from the BCG centre. 
The best-fitting parameters for the galaxy profiles are listed in 
Table |H the corresponding galaxies are indicated in Fig. [5] with 
the same numbering sequence as in Table |4] Snapshots of the 
galaxy- subtracted HST image are shown in Fig. IA.ll 

3.5.3. Mass profiles of cluster galaxies 

In order to reduce the number of parameters to an adequate 
number we need to relate the total mass of cluster galaxies 
to an observable quantity; we will thus assume a scaling rela- 
tion that links the mass of a galaxy to its luminosity. We will 
also assume that the morphological parameters of the underly- 
ing galaxy halo match the values inferred from the galaxy lu- 
minosity; this is a reasonable assumption because of the tight 
correlation between light and mass profile in early-t ype galax- 
ies (see e.g. lGavazzi et al.ll2007l;lKoopmans et al]|2006b . For this 
reason, the galaxy centroid, ellipticity and position angle were 
fixed to the values obtained from the analysis performed with 
SExtractor of the r filter LBC images. However, the parameters 
derived from the LBC data for some of the galaxies close to the 
lensed images or to the BCG appear to be incorrect, because of 
an inefficient de-blending for the lower spatial resolution of LBC 
images. In these cases we adopted the parameters (namely cen- 
troid coordinates, position angle, major/minor axis length and 
magnitude) estimated through SExtractor from the HST F606W 
image. 

We parameterised the mass of cluster galaxies with the 
dual pseudo isothermal elliptic a l mass distribution (hereafter 
dPIE, iKassiola & Kovned 119931: iKneib et al.1 119961) . following 



8 The Tiny Tim software is available 
www. stsci.edu/software/tinytim/tinytim.html . 



at 



the url 



Juli o et al.1 (120071) : iBrainerd et all §996) and several other au- 
thors. The dPIE profile is given by 



Po 



(l + r 2 /r c 2 ore )(l + r 2 /d t )' 



(9) 



The dPIE cut radius can be roughly considered as the half mass 
radiu|J The two-dimensional surface m ass density distributio n 
is obtained by integrating Eq.[9](see also iLimousin et al.ll20Q5h : 



2 

°~ o r cut 



2G(r cu t r core) 



+ R 2 



4 



rout** 2 



(10) 



where <xo is a characteristic central velocity dispersion, which 
cannot be simply rela ted to the measured velocity dispersion 
(lElfasdottir et al.ll2007h . 

The free parameters in the dPIE model are thus the core radius 
r CO re, the truncation radius r cut , and the characteristic velocity 
dispersion <x q. For elliptical galaxies , the Faber- Jackson relation 
(hereafter FJ. lFaber & Jackson! 1 97 6l) predicts that the central ve- 
locity dispersion <xo is proportional to the galaxy luminosity L 
following the scaling relation 



(TO 



I L \ 1/4 



(11) 



where L* is a given arbitrary galaxy luminosity and <Tq is its cor- 
responding dPIE velocity dispersion; our reference magnitude is 
m* = 18.0. We will assume the following scaling relations for 
the characteristic radii of the dPIE profile: 



1/2 



(12) 



here again r* ore and r* ut are the dPIE scale parameters for a L* 
galaxy. The exponent of the scaling relation for r cut determines 
how the mass scales with respect to the light distribution. If 
a = 0.5 the mass-to-light ratio (ML) of the galaxy is constant 
and independent of its lumi nosity, as usually assumed in strong 
and weak lensing analyses (Richard et al.ll2007l: iHoekstra et al.l 



9 The cut radius corresponds to the 3-D half mass radius for a null 
core radius (EKasdottir et al.l l2007l) . 
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12004 . Assuming a - 0.8 one obtains that the m ass-to-light ra- 
tio scales with luminosity accord ing to ML oc L° 3 ( Hal kola et al.l 
2006; Nataraian &JCneib|[l997). In this paper we will assume 
a = 0.5, but in § I3.7l we will compare the results obtained under 
this assumption with those obtained by assuming a = 0.8. We 
did not test alternative scaling relations for the cluster galaxy 
core radius, since we did not identify any reliable constraint on 
the galaxy inner mass distribution. 

The galaxy parameters optimised in the lens modelling are the 
central velocity dispersion and the truncation radius, whereas the 
other dPIE parameters (ellipticity, position and orientation of the 
halo) were tied to the values inferred from the analysis of the 
galaxy light profile. 

The normalisation for the core radius scaling was set to 
r* ore =0.035 arcsec 0.15 kpc, while for the velocity dispersion 
normalisation cr* we adopted the uniform prior [120.0-200.0] 
km s" 1 . Following t he relation between lu minosity and velocity 
dispersion found by Bernardi etaD (l2003h . our reference mag- 
nitude corresponds to a velocity dispersion ^170 km s" 1 for 
an early-type galaxy. The prior cr* e [120.0 - 200.0] km s" 1 
adopts a conservative upper limit; the lower limit was decreased 
to 120.0 km s ~ | to t ake into account the recent findings of 
Nataraia n et al.l (l2009h about the difference in the dPIE veloc- 
ity dispersion between cluster galaxies in the core or in the outer 
regions. 

The truncation radius normalisation determines the steepness of 
the galaxy mass profile. A common and reasonable approach i s 
to fix it to a likely value (see R09 and Limousi n et al . 2009), 
since usually the strong lensing system does not provide enough 
constraints on the mass distribution of the small scale per- 
turbers. However, the mass distribution of early-type galaxies in 
high-density environment is likely to be tr uncated compared to 
field galaxies with an equivalent luminosity (lHalkola et al .120071: 
iLimousin et"aIll2007l:lNataraian et al.ll2002h . For this reason, the 
truncation of cluster galaxy halos could be different from one 
case to another. 

Thanks to the peculiar lensing system of Abell 61 1, an estimate 
of the truncation radius normalisation will be derived in one ded- 
icated tests (see Sect. 13.7b . The perturbing galaxy labelled 1 in 
Fig. [5] was optimised singularly, since the strong constraints on 
its Einstein radius provided by the conjugate systems B.4/B.5 
enable a reliable esti mate of its mass distri bution. Because how- 
ever neither R09 or iNewman et al.l (120091) unlinked this galaxy 
from the cluster member catalogue, we ran a test to verify if and 
how this single perturber can affect the global test results. 

3.6. Modelling of the BCG 

The BCG itself was included in the lens modelling as a dPIE 
halo. However, although the BCG was parameterised with the 
same mass functional as the cluster members, its mass model 
refers only to its baryonic content, while for cluster galaxies 
what is actually modelled is the total mass . As pointed out in the 
work of iMiralda-E scude (1995) (see also ILimousin et al.l l2008: 
iDubinski 1998; Mellie r et al.lll993l) . the dark matter halo of the 
BCG is indistinguishable from the cluster halo in the lens mod- 
elling. 

For this reason a direct comparison of the measured central ve- 
locity dispersion of the BCG to the fitted <To,dpiE would lead to 
incorrect conclusions: the former is determined by the combined 
contribution of the dark and baryonic matter, while the latter is a 
fitting parameter referred to as a "stellar-only" BCG component. 
The dPIE mass parameters - characteristic central velocity dis- 



persion <To,dpiE an d truncation radius r cut - were optimised for the 
BCG in the ranges [100.0-350.0] km s" 1 and [8.0-31.0] arcsec 
(^ [34.6 - 134.2] kpc), respectively: these ranges take into ac- 
count the results obtained through a strong lensing analysis by 
R09 and by N09, to allow a comparison with their results (see 
§0. 

The degeneracy between the central baryonic mass budget and 
the dark matter content with regard to the strong lensing results 
is currently well established. It is still to be discussed whether 
the choice of the BCG mass profile can impact the derived clus- 
ter mass parameters. To assess this point, we ran a dedicated test 
adopting the Sersic profile to parametrise the BCG mass. The 
Sersic mass profile has the same functional form as Eq. [U but 
the parameter H e refers to the mass density at the effective radius 
r e > 

For the modelling of the BCG we adopted the morphological 
parameters derived from the two-dimensional fitting of the BCG 
luminosity profile obtained with Galfit, listed in Table H We op- 
timised the Sersic index within ±10% of the best-fitting value, 
to take into account the measurement errors. After some prelim- 
inary tests, we chose as lower limit to the H e optimisation range 
the value of 2 x 10 8 M©kpc~ 2 to derive a total mass for the BCG 
similar to the mass derived when modelling the BCG as a dPIE 
profile, for a safer comparison between the results obtained in 
both cases (see the total projected mass maps and the surface 
mass density profiles shown in Fig|7]). 




Fig. 7: [Top Panels] Comparison between the projected mass 
maps obtained when modelling the BCG with a dPIE profile (left 
panel, case 1 in TableO or with a Sersic profile (right panel, case 
7 in Table©. See Tableland the text for further details. [Bottom 
Panels] Comparison between the surface mass density profiles 
for the BCG halo [left panel] and the NFW cluster scale halo 
[right panel] derived when modelling the BCG with a dPIE pro- 
file (red solid line, case 1) or with a Sersic profile (green dashed 
line, case 7). 
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Table 5: Median values for the total mass parameters, obtained with the strong lensing analysis. 



Results of the strong lensing analysis of Abell 611 



Case 


x c 


y c 


e 


PA 




C200 


M 200 


cr 




| rmsi 






[arcsec] 


[arcsec] 




[deg] 


[kpc] 




L10 M J 


[km s~ ] 




[arcsec] 




Cluster 


-0.6+0.2 


0.9+0.2 




IJi.OIl/.Z, 


914 ^+9^ fa 


6 78+ - 81 

o./o_ 037 


t.UOIUJ 1 


1 ^00 5+15 7 


0.6 [61] 


0.33 


1 


BCG 

Perturber 1 
Ref. Galaxy 


(0.0) 
(2.33) 


(0.0) 
(-7.85) 






85-5^ 

17 i +6.4 
31 • 1-14.0 






310-9!£o 
176.3^ 

T7Q 1+15.9 

179.1_ 119 








Cluster 


-0.6±0.1 


0.9+0.2 


0.25+0.01 


132.8+0.2 


228.0^-6 


r a +0 56 
6 ' 43 -032 


4.87+0.27 


1307.2+16.4 


0.6 [62] 


0.37 


2 


BCG 

Perturber 1 
Ref. Galaxy 


(0.0) 
(2.33) 


(0.0) 
(-7.85) 






91 1 +84 

i +4.4 
ZJ - A -10.7 

(43.3) 






324 - 7 -33 3 .l 

178.2+6.5 

187. rj^ 








Cluster 


-0.6±0.1 


0.9+0.2 


0.25+0.01 


132. 9+0.2 


o 1 o o +22 7 

218.3+^ 
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Notes. The results are derived assuming an elliptical NFW profile for the cluster-scale lens component and an elliptical dPIE for the galaxy scale 
halos. We report only the values that were optimised in the strong lensing fit; additionally, the coordinates of the cluster galaxies are listed in round 
brackets for the sake of clarity, but they were not optimised. The columns list the coordinates of the halo centroid (expressed in arcseconds with 
respect to the BCG centre: at the cluster redshift and for the assumed cosmology, 1 arcsec=4.329 kpc), the potential ellipticity (here expressed 
as e = (a 2 - b 2 )/(a 2 + b 2 )), the halo position angle, the scale radius r* (that indicates the scale radius r s for the cluster-scale NFW halo and the 
truncation radius r cut for the galaxy-scale dPIE halos), the concentration parameter c 2 oo> the total mass M 2 oo, the characteristic velocity dispersion 
cr for the NFW or the dPIE profiles, expressed in km s -1 , the reduced^ 2 (the number of degrees of freedom is enclosed within square brackets) 
and the root-mean- square computed in the image plane (rmsi). The uncertainties represent the 68% statistical errors, derived assuming a Gaussian 
probability distribution. If the upper and lower errors are significantly asymmetric we report both of them. For a description of the cases listed 
here, see Pag. [12] 



3. 7. Strong lensing results effect of the galaxy parameterisation on the best-fitting results 

for the cluster-scale halo we performed several tests, changing 
We present in this section the results of the parametric analysis the manner of including the cluster galaxies, 
of the strong lensing system in Abell 61 1. In order to assess the 
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Fig. 8: Scatter plots for the NFW parameters r s and C200 obtained for some representative cases. The numbering of the cases corre- 
sponds to the numbering in TableJU The red point distribution marks the 3cr confidence level for the two parameters. The green stars 
indicate the best-fit values; the error-bars refer to the Icr statistical uncertainty. We used the same knot systems as inputs and the 
same priors on the cluster halo mass distribution in all the six cases; the only differences were the modelling and the parameterisa- 
tion of the cluster galaxies and of the BCG. The probability distribution in case 4 appears squeezed because it crushed on the lower 
limit for the BCG characteristic central velocity dispersion. The lens model results for the cases shown here are listed in Table 



We list below the characteristics of the lens modelling for the 
cases presented in Tabj5l 

- Case 1: both the reference truncation radius and the refer- 
ence central velocity dispersion for the cluster-galaxy scal- 
ing relation were optimised in the ranges [4-12] arcsec and 
[120.0 - 200.0] km s"\ respectively. 

In this case, and when not otherwise stated, the mass pa- 
rameters of galaxy n. 1 in Fig J5] were individually optimised 
(the optimization ranges were r cut e [0.5 - 8.0] arcsec and 
cr e [90.0 - 190.0] km s" 1 ) and the exponent a in Eq. [12] 
was set equal to 0.5. 

- Case 2: the reference central velocity dispersion was opti- 
mised, whereas the reference truncation radius was fixed to 
the value of r* ut =10 arcsec. The exponent a in Eq. [T2l was 
assumed to be equal to 0.8 (therefore the Mass to Light ratio 
ML is oc 0.3). 

- Case 3: as in case 2, but we assumed that a = 0.5 (therefore 
ML is constant and independent of the galaxy luminosity). 

- Case 4: a constant ML scaling relation was assumed for all 
cluster galaxies. This is the only case in which the galaxy 
n. 1 in FigUJwas not individually optimised. 

- Case 5: we adopted the same lens model as in case 3, but we 
included as constraint the candidate lens system individuated 



by R09 and by N09, labelled D in FigEJand in TabEOrn the 
Appendix. 

- Case 6: a constant ML scaling relation was adopted for all 
the cluster galaxies, with the exception of those that can be 
tested through the strong lensing analysis, for their direct ef- 
fect on some lensed images (galaxies 1 to 6 in Figj5]), which 
were individually optimised. 

- Case 7: the same lens model as in case 1 was assumed, but 
the BCG was modelled with an elliptical Sersic profile. 

It is not redundant to verify whether the mass distributions 
of the perturbing galaxies differ from the distribution that can 
be derived adopting a scaling relation, since, as shown in 
Table HI these galaxies have significantly different light profiles. 
Moreover, the intrinsic scatter in the FJ scaling relation between 
velocity dispersion and luminosity for early-type galaxies is not 
taken into account in the strong lensing modelling. We stress that 
in all the aforementioned cases, except for case 7, the BCG was 
modelled as a dPIE halo (see Sect. 13.61 for more details on the 
BCG modelling). 

The results of the strong lensing analysis of Abell 611 for the 
aforementioned cases are listed in Table [5] The critical and the 
caustic lines and the images predicted in all tested cases are 
shown in the Fig lB.ll in the Appendix. Evidently depending on 
the different lens modelling the results on the cluster total mass 
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can change by up to 15%, so the parameterisation of the per- 
turber galaxies and of the BCG has a non-negligible impact on 
the strong lensing results; the degeneracy between cluster and 
galaxy mass parameters is shown in Fig 03 In particular, a strong 
degeneracy is found between the fiducial BCG central velocity 
dispersion and the cluster concentration/total mass (see the top 
panels in FigO. 

It is remarkable that the resulting NFW cluster parameters de- 
pend not only on the BCG baryonic mass, but also on the BCG 
mass profile. In case 1 and in case 7 we show the results obtained 
by modelling the BCG with a dPIE or with a Sersic profile, re- 
spectively; the corresponding projected mass maps and surface 
mass density profiles are shown in Fig. [71 The BCG mass de- 
rived in both cases is similar, but the inferred cluster parameters 
are significantly different. In particular, the cluster total mass is 
larger and the NFW concentration parameter lower when mod- 
elling the BCG with a Sersic profile. This result can be explained 
by considering that, in this case, the surface mass density of 
the best-fitting Sersic profile in the BCG central region is larger 
compared with the best-fitting dPIE profile (see the bottom left 
panel in Fig. [7]). The inferred concentration for the cluster scale 
halo in the latter case will thus be higher because of the degen- 
eracy between the cluster halo and the BCG mass with respect 
to the lens total mass, which is the quantity constrained by the 
gravitational lensing. 

However, the estimates of the total mass in the central regions 
(where the strong lensing constraints are found) derived in the 
different cases mutually agree, (see the discussion in Sect.|4]); the 
modelling of the galactic component appears to have a negligible 
impact on the total mass results in the central region. Conversely, 
the differences among several of the test cases arise when the 
cluster mass is extrapolated outside the region that is accessible 
through the strong lensing analysis. 

This is another evidence that the modelling of the BCG stellar 
mass can directly affect the cluster total mass results: in order to 
obtain reliable mass estimates through a parametric strong lens- 
ing analysis, it is crucial to break the degeneracy between the 
different lens components. 

Despite the difference in the best-fitting results, all models return 
probability regions in the parameter space that are marginally 
consistent with each other. As previously found for the X-ray 
analysis, the strong lensing results appears relatively stable 
against possible systematic errors, which in this case are mainly 
constituted by modelling uncertainties. 

As stated before, the perturber galaxy n.l (see Fig. [5]) was in- 
dividually optimised in all the tests that we ran, except in case 
4, due to the tight constraint on its mass distribution provided by 
the conjugated images B.4 and B.5. We found a mild degeneracy 
between the mass of this galaxy, the mass of the stellar content 
of the BCG, and the concentration of the NFW halo, as expected 
given the proximity of this perturber to the BCG (see the bottom 
panels in Fig.©. 

Moreover, in case 6 we optimised singularly the galaxies whose 
mass distribution can be constrained through the strong lens- 
ing effect alone. As it is evident from Fig. [lOl optimising in 
wide ranges the dPIE parameters for these perturber galaxies 
returns, as expected, large uncertainties in the mass parameters. 
However, for some perturbers the parameters inferred through an 
individual optimisation appear to be incongruous with those pre- 
dicted assuming the parameter scaling derived from the Faber- 
Jackson relation. It is reasonable to question whether neglecting 
the scatter associated to the galaxy parameter scaling relations 
in the strong lensing analysis could bias the inferred results. 
It seems likely that in most cases this bias, if any, can be safely 




Fig. 9: Scatter plots of cluster-halo parameters (C200 an d Mot) 
versus parameters of the modelling of either the BCG or the per- 
turber galaxy n.l (hereafter PI, see text for details). Here we 
show the posterior probability distributions obtained in a rep- 
resentative case (case 1 in TabQ. The plotted parameters are: 
1. the NFW concentration parameter C200 versus the BCG fidu- 
cial central velocity dispersion <to,bcg l to P ^ e ft panel]', 2. the 
cluster halo mass M tot versus <to,bcg \f°P right panel]; 3. the 
NFW concentration parameter C200 versus the fiducial central 
velocity dispersion cto,pi of the galaxy PI [bottom left panel]', 
4. the cluster halo mass M tot versus <xo,pi [bottom right panel]. 
The contours show the 68%, 95% and 99% confidence levels. 
The colour code refers to the value of the^ 2 estimator. 



negligible, as it is obvious by comparing case 1 to case 6. 
Nonetheless, in the presence of galaxies that strongly affect the 
strong lensing features, the biasing effect on the final results 
could be more significant (again, this is the case of the perturber 
galaxy labelled n.l in this work, as can be seen by comparing 
case 1 to case 4). 

Regarding the strong lensing systems, all our models pre- 
dict that the tangential arc (system A) originates by the large 
deformation of a single source located on a "naked cusp", 
i .e. a tangential caus t ic extending outside the r adial caustic 
dKormann etal J 119941: iBartelmann & Loebl 119981 iRusin & Mai 
120011) . Naked cusps are a quite rare caustic configuration which 
mainly occurrs in highly elliptical lens systems whose sur- 
face density is only moderately above the critical surface den- 
sity req uired for the production of multiple images ("m arginal 
lenses" iBlandford & Kochanekl 19871: IBa rtelmann 200( 3), with a 
few naked cusps detected so far dLewis et al .1120021: lOguri et all 
2008). Sources placed within the exposed tangential caustic 
form only three magnified images on the same side of the lens: 
this explains why the tangential arc in Abell 611 has no counter- 
image. 

The position and the shape of the observed images are well re- 
produced by all our models, as is obvious from the low scat- 
ter both in the image and in the source plane resulting from the 
best-fit models: all tests returned low values of the reduced^ 2 
0^ 2 ed ^1) and of the image plane root-mean- square (rmsi < 0.5 
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visible [subtracted] in the left [right] panel. 



Fig. 10: Median values of the dPIE truncation radius r cut vs the 
characteristic dPIE central velocity dispersion cr for the galax- 
ies that perturb the strong lensing system of Abell 611. The red 
triangles show the results obtained linking the galaxy parame- 
ters to the same scaling relation in the lens modelling (Case 4 in 
Table [5]) , while the blue stars represent the results obtained op- 
timising the galaxies individually (Case 6 in Table [5]). The grey 
dashed lines show the iso-density contours. The error bars in- 
dicate the 1<t uncertainty, while the cyan boxes indicate the 3<x 
confidence region for the results which refer to Case 4. 



arcsec) (see also the fitting results for a representative case - Case 
1 - listed in Tab©. 

In particular, the lens models of cases 4 and 7 (see Tab© predict 
a central image for the system B in a position almost coinci- 
dent with a radial-oriented feature lying close to the centre of 
the BCG (see Fig Qj] and the predicted central image in Fig jB.li 
this central image was also indentified and interpreted by N09). 
However, since the detection of this feature is not certain, we did 
not include it as an additional constraint. 

Regarding the third system C, our identification of the positive 
parity counter-image, which follows R09, appears reliable ow- 
ing to the very low scatter resulting for this system. Our models 
predict three additional counter-images: two of them should be 
located close to the BCG and to the galaxy n.l, and are pre- 
dicted to be demagnified with respect to the identified image, so 
they are likely not detectable. A negative parity counter-image 
is predicted at about 2-3 arcsec (depending on the different lens 
models) above the observed one. In that area there are no evi- 
dent counter-images, but a faint one, identified with SExtractor 
(S/N 3.9), whose position could be marginally compatible with 
the location of the predicted image. The flux ratio derived when 
modelling the BCG as a dPIE halo does not match the observed 
one well, while the images predicted when modelling the BCG 
with a Sersic profile follow the observed flux ratio more closely 
(see Fig. [T2K because of the difference in the inferred mass for 
the galactic perturbers. This evidence would seem to favour the 
modelling of the BCG with a Sersic profile over the dPIE profile: 
but we aim to highlight that other explanations could be possible 
to the observed mismatch in the flux ratios for the images of the 
third system. For example, the effect of unresolved substructures 
along the line of sight would cause the dim ming of t he negative- 
parity image in cusp systems (lKeetonl2003l) (see also lOguri et al.l 
2008). Another explanation could be the misidentification of the 




Fig. 12: Comparison between the ACS/HST F606W observation 
[left panel], smoothed with a 3 FWHM pixel Gaussian filter for 
a better visualisation of the faint lensed image C.l (marked with 
a white circle), and the reconstructed image systems predicted 
from the lens models in case 1 ([middle panel]: in this case 
the BCG was modelled as a dPIE potential) and case 7 ([right 
panel]: the BCG was parametrised with a Sersic profile - see 
Table 0. The white circles in the middle and right panels mark 
the centre of the predicted counter-images. The black circle in 
the left panel indicates a candidate counter-image of C.l, inden- 
tified with SExtractor. The three images are WCS-aligned; the 
flux scale of the second and the third panel is the same. 



counter-image C.2, which would cause its predicted magnifica- 
tion to be incorrect. 

The candidate system D looks plausible considering the low 
scatter in the observed and predicted positions. The redshift 
estimates of the four knot systems are marginally consistent. 
However, the results of the cluster mass parameters are not sig- 
nificantly changed by including this system. 

4. Comparison between X-ray and strong lensing 
results 

We will compare here the results obtained with the strong 
lensing and the X-ray analyses for the galaxy cluster Abell 611. 
The left panel in Fig. [13] shows the probability distributions of 
the NFW profile parameters obtained with both analyses, while 
the right panel compares the corresponding projected mass 



14 



A. Donnarumma et al.: X-ray and strong lensing analyses of Abell 611 




Fig. 13: [Left panel] Comparison of the NFW parameter probability distribution obtained with our X-ray and strong lensing analyses. 
The X-ray and the strong lensing contours show the 3-cr confidence level (the strong lensing results were marginalised over the other 
fit parameters). The X-ray contours were obtained deriving the X-ray background from the blank- sky data-set, while the strong 
lensing result refers to case 1 and to case 7 in Table [5] For an easier co mparison we oyerplo tted some previous recent results. We 
show in the plot the result s obtained through a st rong lensing analysis by lRichard et a l. ( 2010) (label SL-R09), through a weak and a 
strong lensing analysis by lNewman et al] (120091) (labels WL-N09 and SL-N09, respectively) and in Paper I through a weak lensing 
analysis (label WL-Ro09; the results shown here were obtained by adopting a spherical NFW profile; the shear catalogue in this case 
was derived through the shapelet decomposition technique). The bars show the 1-cr errors. [Top right panel] Projected 2-D mass 
profiles obtained with the X-ray (red solid line) and with the strong lensing analyses when modelling the BCG as a dPIE (green 
dashed line) or a Sersic halo (blue dashed line). The vertical lines mark the spatial range of the strong lensing (black dotted line) 
and X-ray (red solid line) constraints. The X-ray [strong lensing] profile was derived by measuring the enclosed mass in cylinders 
centred on the X-ray emission centroid [BCG]. The distance between the X-ray centroid and the BCG centre is ^ 1 + 1.5 arcsec. 
The mass map derived from the strong lensing analysis was extrapolated up to a radius of 200 arcsec from the mass measurements 
obtained in the central regions. The statistical errors on the strong lensing profiles were determined computing the standard deviation 
from 300 mass map realizations. As a comparison, we report the estimates of the projected total mass within 250 kpc obtained by 
IRichard et al .1120101 through a strong lensing analysis (label SL-R09, marked with a filled black circle) and through an X-ray analysis 
by ISanderson et al.ll2009al (label X-ray-S09, marked with a filled orange square). [Bottom right panel] Ratios between the X-ray 
and the strong lensing projected mass estimates; the strong lensing estimates were obtained by modelling the BCG through a dPIE 
(red solid line) or a Sersic (blue dotted line) profile. 



profiles. For the strong lensing analysis we show the 3cr region, 
marginalised over the other fit parameters, derived in case 1 and 
in case 7 of Table [5j The best-fit values for the NFW profile 
parameters and for the corresponding mass are listed in Table [TJ 
To unfold this issue, higher resolution observations of the 
Sunyaev-Zel'dovich effect measurements towards this cluster 
would play a crucial role, because they would allow us with the 
X-ray data to estimate the projected axis ratio along the line of 
sight. Obviously the best-fit values for the scale radius and the 
concentration parameter derived with the X-ray and the strong 
lensing analyses are consistent within a 2 - cr range; however, 
the results derived from the strong lensing analysis in case 7 
better agree with the X-ray estimates. In particular, the strong 
lensing mass estimates for Abell 611 are lower than the mass 
derived through our X-ray analysis which is lower by up to a 
factor of two. 

This result is somehow surprising, since so far the discrepan- 
cies found between strong lensing and X-ray mass estimates 
have the opposite trend, with the smaller X-ray masses (see 
e.g. IRichard et all 120091: iPeng et al] 120091: iGitti et al] 120071: 
Med ezinski et al] 120071) . However, the mass estimates for 
Abell 611 derived in three representative cases agree very well 
in the region where the strong lensing constraints are found 
(approximately R < 100 kpc), as is obvious from the right panel 
of Fig.[T3j On the other hand, the strong lensing mass is smaller 



than the X-ray one when extrapolating the mass profile in the 
cluster peripheral regions. 

This result confirms that strong lensing provides a robust 
estimate of the total mass of the cluster (DM plus baryonic 
component) in the central region; it moreover suggests that the 
disagreement between the X-ray and the strong lensing masses 
can be related to the incorrect estimate of the total baryonic 
mass budget (for example, the reader can compare the results 
obtained in case 4 and in case 1 in Table [5]) or of its radial 
profile, because the resulting cluster mass distribution will be 
incorrectly steepened (with a higher concentration parameter 
and a smaller scale radius) to account for the missing baryonic 
mass component. 

While this possible source of systematic errors in the strong 
lensing mass estimate was already observed in other studies, the 
dependence of the derived cluster total mass on the assumed 
BCG mass profile is a new evidence in the strong lensing analy- 
sis of real galaxy clusters, which was already observed through 
the a nalysis of simulated galaxy clusters by iMeneghetti et al] 
(2010). A similar but less significant effect is related to the 
scaling of the ML relation, which is assumed to derive the 
cluster galaxy masses, as is evident by comparing case 2 and 
case 3 in Table [5] 

Another explanation of the X-ray/strong lensing mass mismatch 
for Abell 611 could be an oblate-like cluster geometry: as 
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Table 6: Fitting results of the source-plane optimisation. 



Goodness of fit 



System 


x 1 


rms s ["] 


rmsi["] 


A.l 


1.73 


0.067 


0.30 


A.2 


5.56 


0.082 


0.92 


A.3 


2.98 


0.085 


0.40 


A.4 


0.23 


0.024 


0.09 


A.5 


0.15 


0.014 


0.28 


A.6 


1.27 


0.050 


0.12 


A.7 


2.88 


0.066 


0.38 


A.8 


0.98 


0.051 


0.16 


B.l 


2.80 


0.109 


0.29 


B.2 


2.44 


0.100 


0.29 


B.3 


2.37 


0.092 


0.19 


B.4 


2.23 


0.088 


0.22 


B.5 


2.56 


0.083 


0.19 


CI 


0.58 


0.052 


0.10 


Total 


X 1 [d.o.f] 


rms s ["] 


rmsi["] 




28.76 [61] 


0.079 


0.33 



Notes. We show here the goodness-of-fit indicators for a representative 
case (case 1 in Tabled. The columns list the knots system ID, the total 
X 1 for the system, the root-mean- square computed in the source plane 
(rms s ), and in the image plane (rmsi). The last row reports the total x 1 
estimator, referring to a model with 61 degrees of freedom. 

Table 7: Comparison of the mass results obtained with the X-ray 
and the strong lensing analyses. 



Comparison between X-ray and Strong Lensing Results 



[kpc] 



C200 



M200 
[10 14 M o 



X-ray - [Blank sky] 
Strong lensing -/ case 1 ] 
Strong lensing -[case 7] 



350.3 + 79.6 5.18: 



286.9^2 



6.78*$ 



c CQ+0.18 
^'^ -0.31 



0.84 9.3 + 1.4 
4.7 + 0.3 
6.3 + 0.4 



demonstrated in[Gavazzi (2005) (see also Morandi et al. 20 1 Obi: 
aD 120091: Mor andi et al.l l2010at iMeneghetti eTaD 



Corlessetj 

20101) . the orientation of the halo along the plane of the sky 



would cause the lensing mass to under-estimate the true one. 
On the contrary, the X-ray derived masses are the most stable 
against possible pr ojection effect s because of the cluster as- 
phericity (see again Gavazzi 2005). 



5. Comparison with previous studies 

A combined analysis of X-ray data an d of Sunyaev-Zerdovich 
effect measurements was performed by De Fili ppis et al.l (120051) 
to recover the three-dimensional potential distribution for a 
sample of 25 massive clusters. They estimated that the elonga- 
tion of Abell 611 along the line of sight is e\ os ^ 1.05 ± 0.37, 
while its maximum axial ratio is ^ 1.19 ± 0.27, so neither an 
oblate nor a prolate orientation o f Abel l 611 c an be excluded. 
iMorandi. Ettori. & Moscardinil (120071) and iMorandi & Ettoril 
(12007b analysed Chandra X-ray data for a sample of 24 X-ray 
luminous clusters to derive the scaling relations between their 
X-ray properties and the Sunyaev-Zel'dovich effect measure- 



ments (in IMorandi et al.ll2007l) a nd the profiles and the s caling 
properties of their gas entropy (in IMorandi & Et tori 2007). They 
analysed the same Chandra observation of Abell 61 1 as we (see 
Table [T] for more details on the Chandra dataset) by adopting a 
technique similar to the method followed here. Their estimate of 
the total mass of Abell 611 is M 500 = (5.18 ± 2.15) x 10 14 M o at 
the radius - (1.1 ± 0.2) Mpc which agrees very well with 
our X-ray result of M 500 = (6.77 ± 0.84) x 10 14 M o at the radius 
tf 5 oo = (1.2 ±0.1) Mpc. 

A very recent re-analysis of the Chan d ra dat a for 
Abe ll 611 is presented in g ander son et al. (2009a) (see 
also ISanderson & Ponmanl I20T0I: ISanderson et al.l l2009bl) . The 
mass estimate for Abell 61 1 at R^ = (1.3 =b 0.1) Mpc derived in 
their work is M 50 o = (9.23 ± 2.8) x 10 14 M o , which again agrees 
well with our results. 

An analysis of the weak gravit ational l ensi ng effect induced by 
this cluster was performed by iDahlel (120061) through observa- 
tions in the / and V bands obtained with the University of Hawaii 
Telescope and the Nordic Optical Telescope. Assuming an NFW 
profile to model the radial densit y distribution anda concentra- 
tion parameter as predicted by iBullock et al.l (l200ll) . they ob- 
tained a mass value of Mi 80 = (5.21+ 3.47)xl0 14 M o withinfliso . 
The same data were reanalysed by iPedersen & Daniel (120071) . 
who derived an estimate of M 500 = (3.83 ± 2.89) x 10 14 M o . 
A very recent study of the weak lensing effect in Abell 611 is 
presented in Paper I. This work takes advantage of the imaging 
data obtained through the Large Binocular Camera (LBC) 
mounted at the Large Binocular Telescope (LBT), which were 
presented in Sect. Their weak lensing analysis benefits by 
the wide field of view of the LBC (^ 23 arcmin x 25 arcmin) 
and two different methods of e xtracting the shear signal, the 
well known KSB approach dKaiser. Squires. & Broad hursl 
1 19951) and the s hapelet decompo s ition technique, proposed 
mor e recently bv iRefregier) (12003b : iRefregier & Baconl (120031) 
and iMassey & Refregier (2005). Their model-indepent es- 
timates of the total projected mass within a radius of 
1.5 Mpc are M projec ,wL,i = (7.7 ± 3.3) x 10 14 M o and 
M pr0 jec,WL,2 = (8.4 + 3. 8)x 10 14 M o when adopting the KSB or the 
shapelet method, respectively. These results agree very well with 
our X-ray mass estimate, i.e. M pro j ec? x-ray = [9.8 ± 1.4] x 10 14 M© 
at the radius R ^ 1.5 Mpc, while our strong lensing mass 
estimate is lower (see the discussion in the previous section and 
the right panel of Fig. IT3T) . 

They also derived a model-dependent estimate by assuming 
that the total mass profile of the cluster can be parameterised 
through the NFW profile. The best-fitting results were a total 
cluster mass of M 200 = 5.6^[5.9^] x 10 14 M , at the 
radius R 2 oo = [1.57+°-}*] Mpc, and a concentration 



-1.74 



[2.95+} ™] when adopting the KSB 



parameter of C200 
[Shapelets] approach. 
For an easier comparison, the best-fitting values of the cluster- 
scale NFW profile derived in the former case are overplotted in 
Fig- El These results seems to agree better with our X-ray mass 
results as regards the strong len sing estimates. 
Recently. [Richard et al.l (l2010l) (R09) published the results of a 
strong lensing analysis performed on a subsample of 10 massive 
clusters extracted from the Local Cluster Substructure Survey 
(LoCuSS), which included Abell 611. 

Another study dedicated to the analysis of t he matter distr i bution 
in Abell 611 was recently performed by iNewman et al.l (120091) 
(N09). They reconstructed the cluster mass distribution by com- 
bining weak and strong lensing analyses and dynamical data, 
taking advantage of Subaru, HST, and Keck data, respectively. 
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The strong lensing parametric analyses of both R09 and N09 
were performed with the software Lenstool. We overplotted 
their results fo r an NF W profile fit in Fig.[T3l 
iRichard et al j d2010b considered as constraints to the lens 
optimisation both the confirmed (labelled A, B, and C in our 
notation) and the candidate (labelled D) systems (as in the case 
5 in Table 0, but including only one set of conjugated knots for 
the tangential arc (system A in our notation). Their results are 
consistent with our strong lensing findings, which are obtained 
modelling the BCG as a dPIE halo within the 3<x confidence 
range, while a larger mismatch is found with our X-ray results 
and with the strong lensing results in case 7. The discussion in 
§ |H could likely also be applied to a comparison between the 
strong lensing results of R09 and the X-ray and weak lensing 
mass estimates of Abell 611. 

Regarding the results obtained by N09 through the fit of a 
NFW profile to the strong lensing data, their best-fitting values 
show an even higher discrepancy with our mass results: for 
example, they derive a concentration parameter significantly 
higher than the concentration derived from both our X-ray and 
our strong lensing analyses. An extensive comparison of our 
results with the findings of R09 and N09 is beyond the aims of 
this comparative section. One of the differences in the strong 
lensing modelling that could have an impact on the final results 
is the number of adopted constraints, because N09 included 
only the brightest knots of the systems A, B, and D in their 
strong lensing analysis. It could be that using as input all or most 
of the constraints that can be derived from the morphology of 
the lensed images in Abell 611 could decrease the discrepancy 
between their results and ours. Indeed, it is apparent that their 
results follow the degeneracies expected in the strong lensing 
parameter estimates between the cluster halo parameters and 
the central stellar component, which can be reduced by forcing 
the best-fitting lens model to reproduce realistic reconstructions 
of the observed lensed images. Moreover, it is interesting that 
the SL best-fit parameters reported in N09 agree better with 
the results we obtained by linking all the perturber cluster 
galaxies to the same relation (case 4 in Table [5]). This evidence 
reinforces our conclusions about the significant impact that the 
modelling of the galactic component can have in the strong 
lensing analysis of Abell 611, and it could suggest that adopting 
the same scaling relation for the cluster galaxies in this case 
could be inadequate. 

The weak lensing results of N09, instead, agree within the 3<x 
errors with both our strong lensing and X-ray mass estimates. 

6. Conclusions and discussion 

We presented new mass estimates for the relaxed cluster 
Abell 611, derived through the analysis of the cluster X-ray 
emission and of its exceptional strong lensing system. We per- 
formed a non-parametric reconstruction of the gas density and 
temperature profile through the only Chandra observation avail- 
able to date. The cluster total mass was then estimated adopting 
the NFW profile as mass functional, under the assumptions of 
spherical symmetry and of hydrostatic equilibrium. 
Our X-ray estimates of the cluster total mass are M200 = 9.32 ± 
1.39 x 10 14 M Q and M 200 = 11.11 =b 2.06 x 10 14 M Q when deriv- 
ing the X-ray background through the blank-field observations 
or through the source dataset, respectively, thus demonstrating 
the effect of the X-ray background treatment on the cluster mass 
estimate. The inclusion of the metallicity as free parameter in 
the spectral fitting was found to have a negligible impact on the 
cluster mass results. 



We reconstructed the projected mass distribution in the central 
region of Abell 611 through a parametric analysis of its remark- 
able strong lensing system, using the publicly available software 
Lenstool. We used as constraints to the lens optimisation a de- 
tailed system of conjugated lensed images, indentifying 14 sets 
of knots in three different lensed systems. We adopted an ellip- 
tical NFW profile to parametrise the smooth cluster halo; more- 
over we added the small-scale perturbers constituted by the clus- 
ter galaxies by adopting a scaling relation that links their total 
mass to their luminosity. We ran several tests to verify if and 
how the strong lensing cluster mass estimate is dependent on the 
modelling of the galactic component, in particular on the distri- 
bution of the BCG baryonic mass component. 
We found that the extrapolated strong lensing mass for Abell 611 
is 

- moderately sensitive to the mass modelling of the galaxy- 
scale perturbers that determine the peculiar characteristics of 
the strong lensing system in Abell 611; indeed, we obtained 
that the different cluster galaxy parametrisations can induce 
a scatter of up to 10% on the total mass estimate (when 
extrapolating the SL mass up to R200) with respect to the 
mean value (comparing Cases 1 to 6 in Table [5]); 

- more significantly affected by the inferred BCG baryonic 
mass budget (because of the degeneracy between the cluster 
halo mass and the BCG baryonic mass) and by the assumed 
mass profile for the BCG baryonic halo. The deviation of the 
extrapolated total mass with respect to the mean value can 
be as high as 40% when adopting different mass profiles for 
the BCG stellar mass distribution (comparing cases 1 to 7 in 
Table©. 

For example, from the strong lensing analysis we derived for 
Abell 611 the mass estimates of M2oo,dPiE = 4.68 ± 0.31 x 
10 14 M Q , when modelling the BCG with a dPIE model, and 
Sersic — ^-^^_o 23 * wnen assuming the Sersic pro- 

file as the BCG mass functional. The baryonic BCG mass in- 
ferred in both cases is similar, and both estimates are dependent 
on the optimisation limits imposed to the BCG mass parameters. 
Our X-ray and strong lensing mass estimates agree well in the 
central regions, where the total projected mass can actually be 
constrained through the strong lensing analysis, while we find 
a significant difference between the X-ray mass and the extrap- 
olated SL mass in the cluster outer region. This result suggests 
that the observed mismatch between the X-ray and the strong 
lensing mass could be caused by the incorrect estimate of the 
BCG stellar mass content or of the total mass of cluster galaxies 
in the central regions, because gravitational lensing alone cannot 
separate the mass of the different lens components. 
However, even with these caveats in the interpretation of the re- 
sults, the X-ray and strong lensing mass estimates for Abell 611 
in the most opportune cases are consistent within the 3<x error 
range. These findings support the hypothesis that taking into ac- 
count the possible sources of errors and deepening our knowl- 
edge of the current limitations of the X-ray and strong lensing 
analyses, both techniques return robust results. Their reliabil- 
ity can be increased even more by extending the detailed anal- 
ysis of possible systematic biases affecting the mass estimates 
to other appropriate objects, i.e. massive, relaxed galaxy clusters 
like Abell 611. 
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Appendix A: Multiple-image systems 

We list in this appendix the coordinates of the input multiple image systems that we used to constrain the lens model for Abell 611. 
The position of the knots for the image systems corresponding to source A and B are shown in Fig. IA.ll 



Table A.l: Coordinates of the multiple images used as constraints in the lens model optimisation. 



Multiple image systems 



A.l 



A.2 



A.3 



System A 



120.2372 

120.23744 

120.23715 

120.23727 

120.23779 

120.2373 

120.23739 

120.23716 



36.061151 1 
36.060832 2 
36.061185 3 
36.06103 4 
36.060597 5 
36.060982 6 
36.060874 7 
36.061085 8 



120.24078 36. 
120.2392736. 
120.2408736. 
120.24057 36. 
120.23838 36. 
120.2404636. 
120.2399936. 
120.24068 36. 



059618 1 
060178 2 
059567 3 
059648 4 
060427 5 
059611 6 
05979 7 
059646 8 



120.2420636. 
120.2421936. 
120.242 36. 
120.2421436. 
120.2423 36. 
120.24217 36. 
120.24221 36. 
120.2420936 



057857 
056926 
057975 
057605 
056396 
.057434 
05721 
.057782 



B.l 



B.2 



B.3 



System B 



B.4 



B.5 



1 120.2411436.05811 1 

2 120.2411936.058177 2 

3 120.2409 36.058445 3 

4 120.2409936.05842 4 

5 120.24083 36.058612 5 



120.24183 36.055047 1 
120.2419 36.055106 2 
120.24188 36.05483 3 
120.24193 36.054947 4 
120.24197 36.054809 5 



120.2356636.0540891 
120.2356636.0541322 
120.2355436.0540663 
120.2355636.0541114 
120.2354136.0540895 



120.2323 36.061433 1 
120.2323836.061464 2 
120.23237 36.061354 3 
120.2324436.061406 4 
120.2324936.061351 5 



120.236 36.054733 
120.23598 36.054772 
120.2359236.054688 
120.2359 36.054695 
120.2358636.054602 



C.l 



C.2 



System C 



120.24193 36.056063 1 120.2320636.061931 



D.l 



D.2 



D.3 



System D* 



D.4 



1 120.2355136.060743 1 

2 120.23748 36.060537 2 

3 120.2354936.060665 3 

4 120.23561 36.060677 4 



120.2376 36.060501 1 
120.2357 36.060719 
120.2374636.06045 
120.2374236.060507 4 



120.24321 36.053464 1 120.23417 36.055649 



120.24326 36.053421 4 120.23409 36.055685 



Notes. Coordinates are in degrees WCS. The redshif ts of the sources cor responding to System A, System B and System C are za = 0.908 + 0.005, 
z B = 2.06 + 0.02 and z c = 2.59 + 0.01, respectively dRichard et al.ll2010h . 

* System D is not yet confirmed spectroscopically: for this reason, we included it as constraint only in one dedicated test. The redshift estimates 
for the knot groups of this system are listed in Table [5] 
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Fig. A.l: Details of the images used as constraints for the lens optimisation corresponding to the sources A and B (see Table IA.l1 in 
the Appendix). The first three panels on the left show the location of the knots selected along the tangential arc, while the last four 
panels on the right show the conjugate knots identified in the quintuplet s ystem. All panels were selected in the galaxy- subtracted 
image, obtained through the two-dimensional fitting algorithm GALFIT (IPeng et al.ll2002l) . The images were smoothed with a 2 
pixel FWHM Gaussian filter for a better visualisation. 
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Appendix B: Comparison between predicted image configurations 

In this appendix we compare the observed strong lensing system of Abell 61 1 to the lensed images predicted by our best-fit models. 
The top left panel shows the ACS/HST image of Abell 611 (filter F606W), in which the galaxies were subtracted through the software 
Galfit. The last seven panels show the images predicted by our best-fitting lens models, corresponding to the cases listed in Table[5j 




Fig. B. 1 : Comparison between the images predicted by our best-fit models and the observed lensed features. [Top left panel] Galaxy - 
subtracted A CS/HST image of Abell 611. The identified (A,B,C) and candidate (D) lensed systems are marked in green, magenta, 
cyan, and red. [Top middle, top right, middle and bottom panels] The predicted critical (caustic) lines are overlaid in red (blue) on the 
images predicted by the best-fit models listed in Table [5] They refer to a source plane at redshift za = 0.908. The predicted sources 
are drawn in the last seven panels with the same colour as the corresponding image system in the top left panel. The conjugated 
knots used as constraints are shown in Fig. [5] whereas their coordinates are listed in Table |AJ"1 The field dimensions of the panels 
are [35.2 x 35.1] arcsec; all fields are WCS-aligned. The flux scale is the same in the last seven panels. 
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